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The status of lattice calculations of the quark spin, the quark orbital angular momentum, 
the glue angular momentum and glue spin in the nucleon is summarized. The quark spin 
calculation is recently carried out from the anomalous Ward identity with chiral fermions 
and is found to be small mainly due to the large negative anomaly term which is believed 
to be the source of the ‘proton spin crisis’. We also present the first calculation of the 
glue spin at finite nucleon momenta. 
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1. Introduction 

Apportioning the spin of the nucleon among its constituents of quarks and glue 
is one of the most challenging issues in QCD both experimentally and theoretically. 

Since the contribution from the quark spin is found to be small (~25% of the 
total proton spin) from the global analysis of deep inelastic scattering data it is 
expected that the remainder should come from glue spin and the orbital angular 
momenta of quarks and glue. The quark spin contribution from u, d and s has 
been studied on the lattice since 1995 with quenched approximation or with 
heavy dynamical fermions Recently, it has been carried out with light dynamical 
fermions 5.6,7,8 strange quark. We will report the calculation of both the 

connected insertion (Cl) and disconnected insertion (DI) contributions to quark 
spin from u, d, s and c using anomalous Ward identity from the overlap fermion 

As for the quark orbital angular momenta, lattice calculations have been car¬ 
ried out for the connected insertions (Cl) I2.i3.i4,i5,i6,i7,i8,i9^ They are obtained by 
subtracting the quark spin contributions from those of the quark angular momenta. 
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It has been shown that the contributions from u and d quarks almost cancel each 
other. Thus for connected insertion, quark orbital angular momenta turn out to be 
small in the quenched calculation 12,14 nearly zero in dynamical fermion cal¬ 
culations 15.16,17,18,19^ Qj^ other hand, gluon helicity distribution AG{x)/G{x) 
from COMPASS, STAR, HERMES and PHENIX experiments is found to be close 
to zero 20,21,22,23,24^ ^ global fit with the inclusion of the polarized deep inelastic 

scattering (DIS) data from COMPASS and the 2009 data from RHIC gives 
^0.2 

a glue contribution / Ag{x)dx = O.lig Q® to the total proton spin of l/2h with 
Jo.05 

a sizable uncertainty. Most recent analysis of high-statistics 2009 STAR and 
PHENIX data show an evidence of non-zero glue helicity in the proton. For 
Q2 — GeV^, they found the gluon helicity distribution Ag{x,Q'^) positive and 
away from zero in the momentum fraction range 0.05 < x < 0.2. However, the re¬ 
sult presented in ’ has very large uncertainty in the small x-region. Moreover, it is 
argued based on analysis of single-spin asymmetry in unpolarized lepton scattering 
from a transversely polarized nucleon that the glue orbital angular momentum is 
absent Given that DIS experiments and quenched lattice calculation thus far 
reveal that only ^ 25% of the proton spin comes from the quark spin, lattice calcu¬ 
lations of the orbital angular momenta show that the connected insertion (Cl) parts 
have negligible contributions, and gluon helicity from the latest global analysis 
is ^ 40% albeit with large error, there are still missing components in the proton 
spin. In this context, it is dubbed a ‘Dark Spin’ conundrum 31.32 

In this talk, I shall present a complete decomposition of the nucleon spin in 
terms of the quark spin, the quark orbital angular momentum, and the glue angular 
momentum in a quenched lattice calculation. I will then summarize the lattice effort 
in calculating the strange quark spin in dynamical fermions and present a result 
of the total quark spin from a lattice calculation employing the anomalous Ward 
identity and, finally, I will show a preliminary first calculation of the glue spin at 
finite nucleon momenta. 


2. Formalism 

It is shown by X. Ji 33 that there is a gauge-invariant separation of the pro¬ 
ton spin operator into the quark spin, quark orbital angular momentum, and glue 
angular momentum operators 

'^QCD ~ Jq Jg ~ 2^® 

where the quark and glue angular momentum operators are dehned from the sym¬ 
metric energy-momentum tensor 


P = 




( 2 ) 
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with the explicit expression 


Jq — “1“ — f 


- ip + {x X (iD)} Ip 


( 3 ) 


for the quark angular momentum which is the sum of quark spin and orbital angular 
momentum, and each of which is gauge invariant. The glue angular momentum 



X X {E X B) , 


( 4 ) 


is also gauge invariant. However, since it is derived from the symmetric energy- 
momentum tensor in the Belinfante form, it cannot be further divided into the glue 
spin and orbital angular momentum gauge invariantly. 

Since the quark orbital angular momentum and glue angular momentum opera¬ 
tors in Eqs. (3) and (4) depends on the radial vector r, a straight-forward application 
of the lattice calculation is complicated by the periodic condition of the lattice, and 
may lead to wrong results Hence, instead of calculating Jq and Jg directly, we 
shall calculate them from the energy-momentum form factors in the nucleon. 

The Euclidean energy-momentum operators for the quark and glue are 


- 

'mq - 


(- 1 ) 4 ^^/ 

/ 


74 Di +7j Di 


-74 Di - 7 i Di 


V’/, 


n-{^) 


(+*) 


- 2 Xr“*°'^ [Gik Gki + Gik GkiP 

k=l 


( 5 ) 

( 6 ) 


where we use the Pauli-Sakurai representation for the gamma matrices and the 
covariant derivative is the point-split lattice operator involving the gauge link Ug. 
For the gauge field tensor we use the overlap fermion Dirac operator. The 

connection between G^^, and the overlap Dirac operator has been derived 

Tr^ [(Jg,^Dov{x, a;)] = ct Gfj,„{x) + G(a^), (7) 


where Tr^ is the trace over spin, ct = 0.11157 is the proportional constant at the 
continuum limit for the parameter k = 0.19 in the Wilson kernel of the overlap 
operator which is used in this work. The overlap Dirac operator Dov{x,y) is expo¬ 
nentially local and the gauge field Gg^ as defined in Eq. (7) is chirally smoothed so 
that it admits good signals for the glue momentum and angular momentum in the 
lattice calculation 

The form factors for the quark and glue energy-momentum tensor are defined 
as 


{P^s\TlpflgJp,s) = [Tl(-g2)(74p^ -f 7*P4) 


1 X 

-7—T2(-g^)(P4a’iaga + PiCriaQa) -(-'?^)g4<?i 

Zm m 


,(E)i 


u'^-'{p,s). (8) 

Q^9 
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where the normalization conditions for the nucleon spinors are 

s) {p, s) = 1, '^u^^'>{p,s)u^^'^{p,s) = , (9) 

S 


2.1. Sum rules and renormalization 

The momentum and angular momentum fractions of the quark and glue depend 
on the renormalization scale and scheme individually, but their sums do not because 
the total momentum and angular momentum of the nucleon are conserved. We shall 
use the sum rules as the renormalization conditions on the lattice. 

Substituting the energy-momentum tensor matrix elements in Eq. (8) to the ma¬ 
trix elements which define the angular momentum in Eq. (2) and a similar equation 
for the momentum, it is shown that 


L9 = lzLiTm+T2m,,g, ( 10 ) 

{x)g,g = ( 11 ) 


where Z^g is the renormalization constant for the lattice quark/glue operator. {x)q^g 
is the second moment of the unpolarized parton distribution function which is the 
momentum fraction carried by the quark or glue inside a nucleon. The other form 
factor, T 2 ( 0 )g_g, can be interpreted as the anomalous gravitomagnetic moment in 
analogy to the anomalous magnetic moment, ^ 2 ( 0 ) 37,38 _ 


Since momentum is always conserved and the nucleon has a total spin of 


1 

2 ’ 


we write 


the momentum and angular momentum sum rules using Eqs. (1), (10) and (11), as 


{x)q + {x)g = = 1, (12) 

L + J9 = \ [z^ mo) + T2(o)], + [ri(0) + r2(0)], | = ^. (is) 

It is interesting to note that from Eqs. (12) and (13), one obtains that the sum of 
the T 2 ( 0 )’s for the quarks and glue is zero, i.e. 


Z^T^{Q)q + Z^T2{Q)g=0. 


(14) 


We used these sum rules and the raw lattice results to obtain the lattice renormaliza¬ 
tion constants Zg and Zg and then use perturbation to calculate the quark-glue 
mixing and renormalization in order to match to the MS scheme at 2 GeV which 
preserves the sum rules. 


2.2. Results of a lattice calculation with quenched approximation 

Before we present the lattice results, we should point out that the three-point 
functions for quarks which are needed to extract the form factors in Eq. (8) have 
two topologically distinct contributions in the path-integral diagrams: one from 
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r 



t 


(a) (b) 

Fig. 1. Quark line diagrams of the three-point function with current insertion in the Euclidean 
path integral formalism, (a) Connected insertions (Cl), and (b) disconnected insertions (DI). 



connected insertions (Cl) and the other from disconnected insertions (DI) 40,41,42,43 
(See Figs. 1). They arise in different Wick contractions, and it needs to be stressed 
that they are not Feynman diagrams in perturbation theory. In the case of Cl, 
quark/anti-quark fields from the operator are contracted with the quark/anti-quark 
fields of the proton interpolating fields. It represent the valence and the higher Fock 
space contributions from the Z-graphs. In the case of DI, the quark/anti-quark fields 
from the operator contract themselves to form a current loop, which represents the 
vacuum polarization of the disconnected sea quarks. 

It should be pointed out that, although the quarks lines in the loop and the 
nucleon propagator appear to be ‘disconnected’ in Fig 1(b), they are in fact corre¬ 
lated through the gauge background fluctuation. In practice, the uncorrelated part 
of the loop and the proton propagator is subtracted. The disconnected insertion 
(DI) refers to the fact that the quark lines are disconnected. For the nucleon, the 
up and down quarks contribute to both Cl and DI, while the strange and charm 
quarks contribute to the DI only. 

A quenched lattice calculation on has been carried out with 3 valence quark 
masses and extrapolated to the physical pion mass where the numerical details of 
the calculation are given We shall present the results in the following table. 

For the unrenormalized lattice results, we find that [72*(0) + T 2 (0)] (Cl) 
is positive and TKO) negative, so that the total sum including the small 
[T)‘(0) -|- T 2 (0) -I- TKO)] (DI) can be naturally constrained to be zero (See Eq. (14)) 
with the lattice normalization constants — 1.05 and Zg = 1.05 close to unity. 
As discussed in Sec. 2.1, the vanishing of the total r2(0) is the consequence of 
momentum and angular momentum conservation. 

The flavor-singlet which is the quark spin contribution to the nucleon has 
been calculated before on the same lattice We can subtract it from the total quark 
angular momentum fraction 2J to obtain the orbital angular momentum fraction 
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Table 1. Renormalized results in MS scheme at fi = 2 GeV. 



CI(u) 

CI(d) 

Cl(u-l-d) 

DI(u/d) 

DI(s) 

Glue 

(x) 

0.413(38) 

0.150(19) 

0.565(43) 

0.038(7) 

0.024(6) 

0.334(55) 

^2(0) 

0.286(108) 

-0.220(77) 

0.062(21) 

-0.002(2) 

-0.001(3) 

-0.056(51) 

2J 

0.700(123) 

-0.069(79) 

0.628(49) 

0.036(7) 

0.023(7) 

0.278(75) 

9A 

0.91(11) 

-0.30(12) 

0.62(9) 

-0.12(1) 

-0.12(1) 

- 

2L 

-0.21(16) 

0.23(15) 

0.01(10) 

0.16(1) 

0.14(1) 

- 



Fig. 2. Pie charts for the quark spin, quark orbital angular momentum and gluon angular mo¬ 
mentum contributions to the proton spin. The left panel show the quark contributions separately 
for Cl and DI, and the right panel shows the quark contributions for each flavor with Cl and DI 
summed together for u and d quarks. 


2L for the quarks. As we see in Table 1, the orbital angular momentum fractions 
2L for the u and d quarks in the Cl have different signs and they add up to zero, 
i.e. 0.01(10). This is the same pattern which has been seen with dynamical fermion 
configurations with light quarks which was pointed out in Sec. 1. The large 2L 
for the u/d and s quarks in the DI is due to the fact that in the DI is large 
and negative, i.e. —0.12(1) for each of the three flavors. All together, the quark 
orbital angular momentum constitutes a fraction of 0.47(13) of the nucleon spin. 
The majority of it comes from the DI. The quark spin fraction of the nucleon spin 
is 0.25(12) and the glue angular momentum contributes a fraction of 0.28(8). We 
show the quark spin, the quark orbital angular momentum and the glue angular 
momentum in the pie chart in Fig. 2. The left panel shows the combination of u and 
d contributions to the orbital angular momentum from the Cl and DI separately 
while the right panel shows the combined (Cl and DI) contributions to the orbital 
angular momentum from the u and d quarks. 

Since this calculation is based on a quenched approximation which is known to 
contain large uncontrolled systematic errors, it is essential to repeat this calculation 
with dynamical fermions of light quarks and large physical volume. 
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3. Quark spin from anomalous Ward identity 

Attempts have been made to tackle the proton spin decomposition with light dy¬ 
namical fermions configurations. There have been a number of calculations of the 
strange quark spin 5.6,7,8 found the strange quark spin As to be in the range 

from —0.02 to —0.03 which is several times smaller than that from a global fit of 
DIS and semi-inclusive DIS (SIDIS) which gives As « —0.11 The large negative 
contribution from the strange quark is confirmed by a recent analysis ® of the world 
data on inclusive deep inelastic scattering data including COMPASS 2010 proton 
data on the spin asymmetries and the precise JLab CLAS data on the proton and 
deuteron spin structure functions which gives As -I- As = —0.106 ± 0.023 

Such a discrepancy between the global fit of experiments and the lattice calcu¬ 
lation of the quark spin from the axial-vector current has raised a concern that the 
renormalization constant for the flavor-singlet axial-vector current could be substan¬ 
tially different from that of the flavor-octet at the lattice cutoff of ^ 2 GeV. 
The latter is commonly used for the lattice calculations of the flavor-singlet axial- 
vector current for the quark spin. To alleviate this concern, we use the anomalous 
Ward identity (AWI) to calculate the quark spin The anomalous Ward iden¬ 
tity includes a triangle anomaly in the divergence of the flavor-singlet axial-vector 
current 

Nf 

d^Al = 2 ^ mfqfi-f^Qf + iNf2q, (15) 

/=! 

where q is the local topological charge operator and is equal to in 

the continuum. We put this identity between the nucleon states and calculate the 
matrix element on the right-hand side with a momentum transfer q and take the 
Igl —>• 0 limit 

{p's\Af,\ps)Sfj,= \im f^{p',s\2'^mfqfi^5qf + 2iNfq\p,s). (16) 

Lattice theory has finally accommodated vector chiral symmetry, the lack of which 
has hampered the development of chiral fermions on the lattice for many years. It is 
shown that when the lattice massless Dirac operator satisfies the Gingparg-Wilson 
relation + = aDj^D with the overlap fermion being an explicit example 

the modified chiral transformation leaves the action invariant and gives rise to a 
chiral Jacobian factor J = g-2»aTr75(i-5ai3) fermion determinant The 

index theorem shows that this Jacobian factor carries the correct chiral anomaly. 
It is shown further that the local version of the overlap Dirac operator gives the 
topological charge density operator in the continuum i.e. 

Tr75(l - laD,,{x, x)) = ^G“,G“'^‘'(x) + 0{a) (17) 

Therefore, Eq. (15) is exact on the lattice for the overlap fermion which gives the 
correct anomalous Ward identity at the continuum limit. Instead of calculating the 
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matrix element of the axial-vector current derived from the Noether procedure ^0,48^ 
we shall calculate it from the r.h.s. of the AWI in Eq. (15) through the form factors 
defined in Eq. (16). 

In the lattice calculation with the overlap fermion, we note that the renormal¬ 
ization constant of the pseudoscalar density cancels that of the renormalization of 
the quark mass, i.e. Zm Zp = 1 foi the chiral fermion. Also, the topological charge 
density, when calculated with the overlap Dirac operator as in the l.h.s of Eq. (17) is 
renormalized - its integral over the lattice volume is an integer satisfying the Atiya- 
Singer theorem. Thus, when the matrix elements on the right-hand side of Eq. (16) 
are calculated with the overlap fermion and its Dirac operator, the flavor-singlet 
axial-vector current is automatically renormalized on the lattice non-perturbatively 
d la anomalous Ward identity (AWI). 

Besides the fact that AWI admits non-perturbative renormalization on the lat¬ 
tice, the pseudoscalar density in DI and the topological density represent the low- 
frequency and high-frequency parts of the divergence of the axial-vector quark loop 
respectively. It is learned that on the 24^ x 64 lattice, a mere 20 pairs of the overlap 
low eigenmodes would saturate more than 90% of the pseudoscalar loop in configura¬ 
tions with zero modes On the other hand, it is well-known that the contribution 
to the triangle anomaly comes mainly from the cut-off of the regulator. There¬ 
fore, the topological charge density represents the high-frequency contribution of 
the axial-vector loop, albeit in a local form (the overlap operator is exponentially 
local). Since the pseudoscalar density is totally dominated by the low modes, we 
expect that the low-mode averaging (LMA) approach should be adequate for this 
term. To the extent that the signal for the anomaly term is good, we should be able 
to calculate the flavor-singlet qa with the AWI . Both the overlap fermion for the 
quark loop and the overlap operator for the topological charge density are crucial 
in this approach. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 


Q^(GeV2) 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 


Q^(Gev2) 


(a) (b) 

Fig. 3. (a) The charm pseudoscalar and topological density contributions to the proton spin as 

a function of Q^. (b) The same as in (a) for the strange. 
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With the approach described above, we have seen good signals on the 24^ x 64 
lattice with the sea quark mass corresponding to a pion mass at 330 MeV 
We first show the results for the charm quark which contribute only in the DI. 
The pseudoscalar density term and the topological charge density term are plotted 
in Fig. 3(a) as a function of . We see that the pseudoscalar term is large due 
to the large charm mass and positive, while the topological charge term is large 
and negative. When they are added together (black triangles in the figure), it is 
consistent with zero for the whole range of . When extrapolated to = 0, the 
charm gives zero contribution to the proton spin within error due to the cancellation 
between the pseudoscalar term and the topological term. It is shown that the 
leading term in the heavy quark expansion of the quark loop of the pseudoscalar 
density, i.e. 2mP is the topological charge but with a negative 

sign. Thus, one expects that there is no contribution to the quark spin from heavy 
quarks to leading order. It appears that the charm quark is heavy enough so that 
the 0{l/'m?) correction is small. We take this as a cross check of the validity of our 
numerical estimate of the DI calculation of the quark loop as well as the anomaly 
contribution. 

The contributions from the strange are also calculated and shown in Fig. 3(b). 
The 2mP contribution is slightly smaller than that of 2q and results in a net small 
negative value for the sum of 2mp and 2q at finite . After a dipole fit, we obtain 
As = —0.026(5) at = 330 MeV. Here, As denotes the contributions for both s 
and s. Alt and Ad are similarly defined in the following. 


mq=500 MeV light iso-vector (m,j=330 MeV) 




(a) (b) 

Fig. 4. (a) The quark spin of the proton-like baryon with mq ~ 500 MeV from both the axial 

vector current and the pseudoscalar term through AWL In this case, the DI contribution of 2mP 
is canceled by the topological charge term, (b) The same as in (a) for light quarks at the unitary 
point for isovector which involves only Cl. 


Since this As is quite a bit smaller than the experimental value, we explore the 
possible finite volume effect and the fact that the induced pseudoscalar form factor 
hA{q^) has been neglected in the extrapolation which does not contribute at 
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the = 0 limit as in Eq. (16), but has a contribution at finite We shall 

check this in the connected insertion (Cl) calculation. As can be seen in Eig. 4 for 
niq ^ 500 MeV, both Au and Ad in Cl calculated from the axial-vector current 
and renormalized with Za from the isovector Ward identity are well reproduced 
through the extrapolation of 2mP with a dipole form. Whereas, in the case 
of light quarks at the unitary point, g\ = 1.13(2) from the axial-vector current 
is 1.8(1) times larger than 0.62(4) from the dipole extrapolation of 2mP. This is 
most likely due to the ignorance of the induced pseudoscalar form factor hA{q^) as 
well as the finite volume effect at small which is well known to plague the 
extrapolation of the nucleon magnetic form factor. 


light (m^=330 MeV) 



Fig. 5. The combined pseudoscalar contribution from both the connected insertion (Cl) and DI 
in the plot), along with the overall quark spin from both pseudoscalar and topological 
charge [qa d)- The plot corresponds to the unitary point with m,r = 330 MeV. 


At the unitary point, when the valence u/d mass matches that of the light sea. 
Fig. 3 shows the quark spin contribution from the combined pseudoscalar terms 
2mPu/d of the Cl and DI with a dipole extrapolation. Also plotted are the overall 
quark spin Au/Ad by including the topological charge contribution. In this case, 
we obtain Au -\- Ad = 0.19(3) and Au — Ad = 0.62(4) at = 0 from a dipole 
extrapolation in . As we discussed above, the fact that g\ from the axial current 
is 1.8(1) times larger than that of Am — Ad through the Ward identity approach 
is most likely due to the neglect of the induced pseudoscalar form factor hA{q^) 
and the finite volume effect in the extrapolation. We apply this 1.8(1) factor as 
an estimate to correct the present AWI approach and obtain Au -I- Ad = 0.35(6), 
As = —0.05(1). Thus the total estimated spin AE = 0.30(6) at the unitary point 
is consistent with the present experimental results which are between 0.2 and 0.3. 
We expect that, at lighter quark masses, AS will be smaller. 
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The above results are from the 24^ x 64 lattice at = 330 MeV with 200 
configurations. The nucleon propagator in the DI has been calculated with the 
smeared-grid noise source with time dilution which covers all time slices in order to 
have reasonable statistics for the DI. 


4. Glue spin 

It has been pointed out that decomposing glue angular momentum into glue spin 
and orbital angular momentum is only feasible in a specific gauge Making contact 
with the parton picture, a spin sum rule involving quark and glue spins and orbital 
angular momenta is derived in the light-cone gauge (i.e. A~^ = 0) with nucleon in 
the infinite momentum frame The longitudinal glue spin content is 

'S'g = (p,s| J d^xTr{ExAf\p,s)/{p,s\p,s), (18) 

where the nucleon state is in the infinite momentum frame and the gauge potential 
and the gauge field are in the light-cone gauge. Similarly, a gauge-invariant glue 
helicity distribution is defined with the light-cone correlation function 

^ J ^e-“^-"«-(P5|F+“(r)>C“''(r,0)F„+(0)|P5) , (19) 

where = (1/2) is in the adjoint representation with A'^ = T'^A'^, so 

is the light-cone link £(^”,0) = Pexp[—ig fj ,0±) dr]~] . 

Since lattice QCD is formulated in Euclidean time, it is not equipped to address 
the light-cone gauge or the light-cone coordinates and; as such, one is not able to 
calculate AG as defined in Eqs. (18) and (19) on the lattice directly. 

On the other hand, a gauge-invariant decomposing of the proton spin has been 
formulated and examined in various contexts 58,59,60,61 _ jg leased on the 

canonical energy momentum tensor, instead of that in the symmetric Belinfonte 
form. The glue spin operator is 

Sg = E^x Al^y, ( 20 ) 

where A^ phys is the physical component of the gauge field Ap which is decomposed 
into Apphys and a pure gauge part as in QED, 

Ap Appfiyg -\- App^j.Q. (^ 1 ) 


They transform homogeneously and inhomogeneously with respect to gauge trans¬ 
formation respectively, 

Apphys -^fiphys ~ 9-^p.physg 

Appure -^ypure ~ gAppureQ ~ ~9^p9 ) ( 22 ) 

where g is the gauge transformation matrix and go is the coupling constant. In oder 
to have a unique solution, conditions are set as follows: the pure gauge part does 
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not give rise to a field tensor by itself and Aphys y, satisfies the non-Abelian Coulomb 
gauge condition 


Fyu pure — ^yAjj pure ^i^Ap pure '^doiAppure: A^, pure] — 0 

^iAiphyS — diA^phyg A^p^yg] — 0 . (^ 3 ) 


This is analogous to the the situation in QED where the photon spin and orbital 
angular momentum can be defined 62,63,64,65 fj-Qjjj canonical energy-momentum 
tensor 


Sa 


E± X A± 


(24) 


La = 



X X V) A^ d^x, 


(25) 


where _L denotes the transverse part. Since they are defined in terms of the trans¬ 
verse parts, they are gauge invariant. However, this gauge invariant definition breaks 
Lorentz invariance. Nevertheless, it is shown that the ‘spin’ and ‘orbital’ angular 
momentum so defined are conserved for a free field Eurthermore, they are ob¬ 
servables and can be measured in experiments through interaction with matter. In 
1936, Beth had observed one component of the spin angular momentum of light 
by measuring the tongue on a birefringent plate exerted by a circularly polarized 
light. Also, it is shown that the orbital angular momentum of a paraxial laser 
beam can be measured. Even though gauge invariance is preserved in this canonical 
formulation, the spin and orbital AM operators are not boost invariant. Since the 
experiments are conducted in the lab, the formulation is adequate for this single 
reference frame. 

After integrating the longitudinal momentum x, the light-cone operator for the 
matrix element has the following expression for the glue helicity 


Hg = 


E;“(0) X (A“(0) - —(VA+’'’)£'’“(r,0)) 


(26) 


It is recently shown that when boosting the glue spin density operator Sg in 
Eq.(20) to the infinite momentum frame (IMF), the second term in the parentheses 
on the right side of Eq. (26) is Apure- Thus Hg is the glue spin density operator Sg 
in the IMF along the direction of the moving frame. In other words, the longitudinal 
glue spin operator turns into the helicity operator in the IMF. 

To carry out a lattice calculation of the matrix element of the glue spin operator, 
it is realized that Apppys is related to that fixed in the Coulomb gauge, i.e. 
Apphya = g~^Augc where A^ is the gauge potential fixed to the Coulomb gauge and 
5 c is the gauge transformation that fixes the Coulomb gauge. Since Sg is traced over 
color, the spin operator is then 


Sg = / d'^xTr[gcEg^ ^ x Ac) = / d'^xTr{Ec x Ac) 


(27) 
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Fig. 6. The results of glue spin Sq in longitudinally polarized proton with longitudinal momenta 
at 0, 460 MeV and 930 MeV. The quark masses in the nucleon propagator correspond to m,r = 380 
and 640 MeV. 


where Ec is the electric field in the Coulomb gauge. Although it is gauge invariant 
since both E and Aphys transform homogeneouly, it is frame dependent and thus 
depends on the proton momentum. Its IMF value corresponds to AG which is mea¬ 
surable experimentally from high energy proton-proton scattering. The important 
outcome of the derivation is that glue spin content is amenable to lattice QCD cal¬ 
culation. To the extent that it can be calculated at large enough momentum frame 
of the proton with enough precision, it can be compared to the experimental glue 
helicity AG. 

The first attempt to calculate Sq on the lattice has been carried out on the same 
set of 2 -I- 1 flavor dynamical domain-wall configurations on the 24^ x 64 lattice 
with the sea pion mass at 330 MeV The electric field E is constructed from 
the overlap Dirac operator defined in Eq. (7). The gauge potential A is obtained 
from the unsmeared gauge link. We obtained results for the longitudinal nucleon 
momenta Pz = n(27r/La) with n = 0,1, 2 which correspond to 0, 460 MeV and 920 
MeV and for the case of quark masses in the nucleon propagator which correspond 
to m-TT = 380 MeV and 640 MeV. The unrenormalized results are presented in Fig. 6. 
We see that the preliminary results in Fig. 6 are quite noisy and, as a result, one 
cannot discern the Pz behavior. The signal can be improved by smearing the link, 
but it is a challenge to reach large Pz on the lattice. Since one needs Pz to be less 
than the cutoff, i.e. PzO, ^ 1 to avoid large discretization error, this will require a 
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large lattice size L so that m.„La > 6 for the nucleon. We note that the question 
how large a Pz is needed to have the quasi PDFs coincide with the PDFs has been 
studied in a spectator diquark model . It is found that it is necessary to have Pz 
as large as 4 GeV for the quasi PDFs to be a good approximation of the PDFs. 

5. Summary 

We have reported the current lattice efforts in calculating the quark spin, quark 
orbital angular momentum, glue angular momentum and glue spin in the nucleon. 
A complete decomposition of the proton momentum and spin into its quark and glue 
components is given in a quenched approximation. In this case, the glue angular 
momentum is not further divided into spin and orbital angular momentum parts. 
The quark spin calculation is recently carried out from the anomalous Ward identity 
with chiral fermions and is found to be small mainly due to the large negative 
anomaly term which is believed to be the culprit of the ‘proton spin crisis’. An 
exploratory lattice calculation of Sq in the non-Abelian Coulomb gauge is carried 
out which has large errors and the nucleon momentum is limited to ^ 1 GeV. 
The signal of the glue spin Sq can be improved with smearing, but the major 
challenge is to have a lattice with fine enough lattice spacing to accommodate large 
momentum states and show that the infinite momentum extrapolation can be made 
under control. 
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